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Using an exactly solvable model, we study low-energy properties of a one-dimensional spinless 
electron fluid contained in a quantum-mechanically moving wire located in a static magnetic field. 
The phonon and electric current are coupled via Lorentz force and the eigenmodes are described by 
two independent boson fluids. At low energies, the two boson modes are charged while one of them 
has excitation gap due to back-reaction of the Lorentz force. The theory is illustrated by evaluating 
optical absorption spectra. Our results are exact and show a non-perturbative regime of electron 
transport. 
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Electronic and mechanical properties of nano- 
electromechanical systems (NEMS] are recently of fun- 
damental and technical interest yj. Perpetual trend 
of miniaturization of electronic devices has already led 
NEMS smaller down to molecular scale, e.g., nanotube 
electromechanical systems which are likely the bases 
of new novel devices in NEMS application^. These sys- 
tems are close to the regime of quantum-mechanical op- 
eration which is the ultimate limit of nano electromechan- 
ical devices. For the design of quantum mechanically op- 
crating NEMS devices, it is crucial to understand first the 
quantum effects of mechanical motion coupled to electric 
current. Here we ask ourselves; Are the electric currents 
influenced by the quantum nature of mechanical motion? 
If so, what is the signature of the quantum fluctuation? 

To answer these, we study one-dimensional spinless 
electron fluids coupled to their quantum-mechanical mo- 
tion in a static magnetic field. This model is exactly 
solvable as will be shown here. In this way, we could 
describe both the electrical and mechanical degrees of 
freedom in a purely quantum mechanical manner on an 
equal footing. The harmonic-fluid approach used for 
electron fluids in a static wire, which are referred to 
as (Tomonaga-)Luttinger liquid theory 0,0,0, could be 
successfully applied to describe the electric currents cou- 
pled to mechanical motions through Lorentz force. This 
magneto-electro-mechanical coupling is essential in so- 
called electromotive-force technique which is popularly 
used in NEMS-based force sensors0Q- 

In this paper, we show that electric currents in nano- 
electromechanical wires are strongly influenced by the 
quantum fluctuations of mechanical motion through 
the magneto-electro-mechanical coupling. This coupling 
causes interesting and peculiar phenomena because of the 
feedback reaction between electric and mechanical move- 
ment. As we shall see, an important frequency scale re- 
lated to a gap structure arises due to this feedback reac- 
tion. The feature we found could be observed generally in 
other vibrating one-dimensional electronic systems such 



as atomic wires 9] and quasi-one-dimensional metals [lfjj. 

We model the electrons with an interacting spinless 
Fermion fluid with a periodic boundary condition ^(x + 
L) = $>(x). These Fermions are contained in a wire of 
mass density p m and tension T, and a magnetic field B is 
applied perpendicular to x-axis. The displacement u(x) 
is in the direction of x x B, and assumed small (d x u <C 1) 
and periodic (u(x + L) = u(x)). The electrons and the 
wire are assumed strongly tightened so any force exerted 
on the electrons in u-direction is delivered to the wire 
immediately. The Hamiltonian reads (h = 1, e < 0) 
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Here n u is the conjugate momentum of u[x); 
[u(x), n u (x')] — i5(x — x 1 ). A x = Bu(x) + &b/L is the 
gauge potential for the applied magnetic field where $s 
is a constant magnetic flux enclosed by the circle of ra- 
dius L/2-k. For the sake of convenience, let us redefine 
u(x) + <&b/BL — > u{x) so that A x = Bu(x). 

Following the approach of Haldane|6j, we use a Bo- 
son field 9 defined from the Fermion density p(x) = 
(k F + d x 9)/n and its conjugate field <f) ( [4'{ x )^{ x ')] = 
i^nsgx\{x — x')). The boson fields describe the long wave- 
length fluctuations of the Fermion field w(x). It proves 
convenient to define new Fermion field ip' = exp(«A)\E r T j 
where A(x) = eB J u(x')dx' . In terms of 8,(j> and A, the 
new Fermion field vb^ is written as 



vb\x) 



exp[in[6(x) + k F x}} exp[i((f>(x) + A(x))} (2) 
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The boundary conditions for the Boson fields are 9(x + 
L) = 0{x) + ttN, <j>(x +L) = <j>(x) + ttJ, and A (a; + L) = 
A(x) + eBuoL, where uq — ^ J Q dxu(x), and the sum of 



the topological numbers N and J is an even integer. Note 
that the new Fermion field now has a twisted boundary 
condition ip(x + L) = cyLp(—ieBuoL)i^(x). 

The effective Hamiltonian describing long-distance 
physics is given by taking terms of n = ±1 in Eq.J2J 
with a proper normalization constant; 



H = J dx{g(d x <j> + eBu) 



H u . (3) 



Here the dimensionless constant g is governed by 
electron-electron interaction strength and v — vp/g is 
the electronic wave velocity, where vp is the Fermi veloc- 
ity for the non-interacting case. For short-ranged elec- 
tron interactions V(x — x') = VS(x — x'), g is given by 
1 / \J\ + V/nvp. In terms of the boson fields, the electri- 
cal current density is given by 



J{x) — —evpd x 4>{x) — vpe 2 Bu(x). 



(4) 



Note that the term of u 2 in Eq.@ and H u describe 



phonons of finite mass at small momenta; 

dxu 2 (x) + H u = ^2u s {q)a\a q , 



^e 2 B 2 
2tt 



Here 



v 



B ' 



(5) 



(6) 



where v s = \/T/p m is the bare sound velocity and 
tug = ^—B. The first term in Eq.JSJ) essentially gives 
mass to phonon and originates from back-reaction pro- 
cess in the system. Note that when the wire moves, 
an induced electric field £ Bm i = —Bdtu changes the 
momentum of electrons by J dte£ cm [ ~ —eBu, which 
causes the change in electric current —e 2 Bupo/m where 
Po = kp/ir is the mean Fermion density. The Lorentz 
force —vpe 2 B 2 u/n exerted on the induced electric cur- 
rent again acts on the movement of the wire, which we 
call here the back-reaction process. 

Now let us use the following Fourier-transformed field 
operators; 



9(x) 
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u(x) 
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Here (#o,J) and (0o,N) are pairs of conjugate action-angle 
variables [6| and [mo,7I"mo] = i- 
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FIG. 1: The dispersion relation of the eigenmodes u±(q) 
(solid lines). u> is in unit of ujb and q is in unit of lub/v. 
The dotted lines denote the boson mode for electron ui — vq, 
phonon ui = v 3 q, and the massive phonon ui s (q). We chose 
v s = 0.3v in this figure. 



Using these operators, the Hamiltonian in Eq.© is 
written as 

H = }^v\q\blbg +}^Q s (q)ala q + H ZM 
q i 



Here Hzm describes the zero-mode; H 



ZAI 



2L *> J 



^f^) 2 + 2^fr + 5^ n2 ' according to which the center- 
of-mass of the wire executes a simple harmonic oscillation 
with resonance frequency ujb- 

The Hamiltonian in Eq.JS|)is now diagonalized using 
new Boson operators 7 9 ,±, which are linear combinations 
of the four boson operators b q , b*_ q , a q , and a}_ q ; 

H = J2 w ±(9)7g,±7g,± + h zm (9) 
q,± 

where the dispersion relations for the two eigenmodes are 
UJ ± (q) = (10) 



(v 2 + v 2 )q 2 +uj 2 b ± ^/{{v 2 - v 2 )q 2 - uj 2 b ) 2 + Auj 2 B v 2 q 2 



Now the ground state is the vacuum of new bosons, 
7 9 ±|0>=0. 

For detailed discussion, let us assume v > v s which is 
usually the case. In Fig. 1, we plot dispersion relations 
for the two modes. As q increases or B decreases, the 
system goes to a weak coupling regime where the '+' 
mode becomes electronic (7 ? + ~ b q , tu+(q) ~ vq) and 
the '— ' mode becomes phonemic (7,- ~ a q , oj-{q) ~ 
u> s (q)). In the opposite limit where q is small and B ^ 
0, the system is in a strong coupling regime where the 



two uncoupled modes undergo level crossing as shown in 
Fig. 1 and eventually completely exchange their identity 
in q — y limit, i.e., 7 g + — > a q and 7 g _ — » 6,j. It is 
instructive to note that electric current can flow through 
the '+' modes even in q — > limit, where the '+' modes 
become phononic. This is because of the 'diamagnetic' 
contribution of electric current —vfs 2 Bu(x) in Eq.Q. 

In the strong coupling regime, the dispersion relation of 
the '+' mode is in the same form with plasmon dispersion 
relation, u)+(q) ~ \J v 2 q 2 + lo 2 b (q <C ojb/v). In contrast 
to the plasmon^]] j however, the frequency gap lob origi- 
nates from the back-reaction which gives mass to phonon 
at long- wavelength limit, Co s (q) — > ujg. The ' mode is 
gapless LO-(q) ~ v s vq 2 /lob (q *C lob/v s ). This mode is 
also charged at small g, where electronic and phononic 
excitation are entangled with each other. These excita- 
tions at long wavelengths are non-perturbative and can 
not be captured by perturbative expansions in B. 

Now let us consider experiments where the mechanical 
quantum fluctuations can be observed. In Fig. 2, we 
depict a system we consider. An external electrical field 
V e (t) I d is applied along the wire in < x < d by nearby 
external metallic gates with oscillating voltage difference 
V e (t) — V e (uj)e iu ' t . If the current I(to) is measured far 
away from the regime < x < d, the current will not be 
sensitive to the position of the contacts so can be given 
by its averaged value I(u>) ~ ^ dxJ(x,uj). From the 
linear response theory, one gets 



a(ui) 



in 



V e (uS) iujd 2 



dx / dx'U(x,x',Lu),(n) 



where U(x,x'co) = i J °° ds < [J(x,0),J(x',-s)] > e tuJS 
is the current-current correlation function. 

This form is equivalent to the conductance formula 
G(u>) used by Kane and Fisher^] but caution is needed 
here. In our case, additional electromotive force V em { (u>) 
is induced along the wire when B =/= 0. In usual two- 
terminal conductance measurements, since the measured 
voltage V{lo) is given by V e mi(w) + V e (ui), the two- 
terminal conductance is related to our a(u>) as I(lo) — 
G(w)(V e (w) + VUfM) = o(«)V e (w). 

The real part of a(ui) should be understood as an 
optical absorption spectrum. After averaging over a 
time period 2it/lo, the power P(u>) absorbed from the 
external electric field to the wire is given by P(u>) = 
±Retfdx(V e *(u J )/d)J(x,u J ) = ±Rea(u)\V e (u)\ 2 . 

To evaluate conveniently the current correlation func- 
tion, we switch to a Lagrangian formalism following Kane 
and Fisher 0| . The Hamiltonian in Eq.© is equivalent 
to the following Euclidean action, 



+ ZJ drdx (^d x u) 2 + ^{d T uf ) (12! 
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FIG. 2: Schematic figure for the experimental set-up in con- 
sideration. The wire is on an insulating layer above two metal- 
lic gates denoted by M. 
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where J(x, r) = ied T 9(x, r) / n. Evaluation of the current 
correlation function Il(x,x',iuj n ) at Matsubara frequen- 
cies io„ = 2im/(3 can be easily done using the above 
Euclidean action; 

Il(x, x', iui n ) = / dre- l ^ T < T T J(x, t)J(x',0) > 
Jo 

/ UJ 2 

dqe iq(x -* > — — ^ - — .(13) 
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Now, by performing analytic continuation H(x,x',lu) — 
H(x, x', iu> n —y uo + + ), we get Rea(w) from Ea. l|ll|l : 

Rea(w) 

_ I I f°° a sin 2 (grf/2) lo 2 - (qv s ) 2 

^ VFM ]_^ aq (qd/2) 2 w 2 -(iO 2 + (q) + L0 2 _(q))/2 

(14) 
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The result of the above integration is written in terms 
of dimensionless parameters Co = \u>\d/v s , B = Llob/v s , 
and r\ = v s /v; Kea(to) — g%- J2± A±(u>, B, rj). (Here we 
revived the Planck constant.) The dimensionless function 
A± are 



A+ = ± 



(15) 



sui 2 (k ± /2) i](Co 2 - k±)/k± 
( K ±/2) 2 J(l- 7 f) 2 Cb 2 +4B 2 ' 



where k± = J ^(1 + ri 2 )Co 2 =p Co\J (1 - r] 2 ) 2 Cb 2 + AB 2 ^j /2. 

In Fig. 3, we plot the field absorption power Rea(w) 
using the exact formula we obtained. Since usually v s <C 
v (for metallic carbon nanotubes, v s ~ 15 km/s <C v 
and v ~ 10 5 km/s ), let us assume i)« 1. In this case, 
the induced electric current is strongly suppressed when 
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with the back-reaction 
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FIG. 3: The optical absorption spectrum for v s /v = 0.1 and 
ujbcI/v = 10 -4 . The inset shows the details at low frequencies. 
In the low frequency regime (inset), the maximum arises near 

id W (v s /v)lUB- 



\lj\ < ujb- There are three interesting frequency ranges, 
(i) lj <C (v s /v)ojb- Here only '-' modes arc excited which 
are mainly electronic density fluctuations. In this regime, 

using A_ » |\/¥ < 1 and A + = 0, we get 



Kea(oj) ~ g 




(16) 



It is interesting to note the anomalous frequency and 
magnetic-field dependence Rea(w) oc y/uj / B which 
can not be obtained within a perturbation analysis, 
(ii) (v s /v)ujb <C u> < lob- Here also only the '— ' 
modes are excited while their nature becomes partially 
phonon-like as the frequency increases. In this regime, 
A_ w 4i]B 2 sm 2 (Cj/2) and A+ = 0. In particular for 
lj < 1 (i.e., lj < v s /L), we get 



UJB 



P> 



.fin 



B = 0.88 x 10 _2 S[T]^ 



' vf[iti/s] 
p m [kg/m 



Hz(19) 



where we recovered h. As clear in its form, electron- 
electron interactions do not affect this characteristic fre- 
quency. For single- walled carbon nanotubes, the linear 
mass density is estimated to be p m ~ 4 x 1CU 13 kg/m 
when its radius is ~ 4 nm . Using a typical Fermi ve- 
locity vf — 10 5 to/s, we estimate uj b ~ 4.4 B[T] MHz, 
which is in the experimentally accessible range. 

While we consider electro-mechanical coupling via 
Lorentz force, capacitive coupling to external gates |2lll3| 
might be of interest in other contexts. In this case, the 
coupled charge density and phonon displacement is de- 
scribed by ~ J dxu(x)d x Q . We assumed no external 
source of mechanical energy dissipation of the wire, which 
is another subtle issue of electromechanical system [hi] 
and important to clarify the detailed operating condi- 
tions of quantum NEMS. 

In summary, we outlined a theory for electrons in one- 
dimensional wires in quantum mechanical motion. Quan- 
tum mechanical displacement coupled to electric cur- 
rent and the back reaction from electrons give rise to 
a current-carrying boson mode with an energy gap. The 
induced electric current by oscillating electric field might 
reveal the elementary excitations in the system through 
its frequency dependence. We presented here a non- 
perturbative regime where the induced electric current 
has non-trivial frequency dependence. 
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(iii) Finally when ojb < oj, the '+' modes play a signifi- 
cant role and A + has singularity near u> — ljb ■ Near this 
singular point, in spite of the phononic nature of the '+' 
modes, significant electric current may flow through its 
'diamagnetic' contribution; 



Rea(a>) ~ g- 



(18) 



The singularity of a{u>) originates from the density of 
states of the massive boson mode 0J+(q) near oj = ljb- 
Similar to Landauer conductance quantization, a quan- 
tization rule Reev((jj) « ge 2 /h appears in weak coupling 
regime ljb -C lj < v/d. 

Important to all the physical consequences discussed 
in this work is a characteristic frequency ljb associated 
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